A double-ended queue with catastrophes and repairs, 
and a jump- diffusion approximation* 



Antonio Di Crescenzo*^^-' 
Virginia Giorno*^^) 
Balasubramanian Krishna Kumar*^^^ 
Amelia G. Nobile^^) 
(1) Dipartimento di Matematica e Informatica, Universita di Salerno, 
Via Ponte don Melillo, 1-84084 Fisciano (SA), Italy 
email: {adicrescenzo, giorno, nobile}@unisa.it 
(2) Department of Mathematics, Anna University, Chennai 600 025, India 
email: drbkkumar@hotmail.com 



Abstract 

Consider a system performing a continuous-time random walk on the integers, subject 
to catastrophes occurring at constant rate, and followed by exponentially-distributed repair 
times. After any repair the system starts anew from state zero. We study both the transient 
and steady-state probability laws of the stochastic process that describes the state of the 
system. We then derive a heavy-traffic approximation to the model that yields a jump- 
diffusion process. The latter is equivalent to a Wiener process subject to randomly occurring 
jumps, whose probability law is obtained. The goodness of the approximation is finally 
discussed. 
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1 Introduction 



The description of numerous types of systems subject to random evolution can be rendered in 
a more effective way by including the eventual occurrence of catastrophes. This is well known 
in the contexts of population dynamics (see, for instance, the contributions by Brockwell 1985, 
Swift 2000 and 2001, and Gani and Swift, 2006 and 2007). In this research area, Chao and 
Zheng (2003) considered an immigration birth and death process with total catastrophes and 
studied its transient as well as equilibrium behavior. Other contributions are due to Economou 
and Fakinos (2003) and (2008), that extend some previous results on birth-death processes with 
catastrophes to the more general cases of continuous-time Markov chains with catastrophes re- 
sulting from a point process, and of the non-homogeneous Poisson process with total or binomial 
catastrophes. Various results on a class of stochastic models for systems subject to random reg- 
ulation via disasters are presented in Stirzaker (2006) and (2007). Also the analysis of queueing 
systems or birth-death processes subject to catastrophes has attracted the attention of several 
investigators (see, for instance, Krishna Kumar and Arivudainambi 2000; Krishna Kumar and 
Pavai Madheswari 2003; Di Crescenzo et al. 2008). A suitable approach adopted by Krinik et 
al. (2005) to determine the transient probability functions of some classical queueing systems 
with catastrophes uses dual processes, randomization and lattice path combinatorics, whereas a 
case of finite birth-death processes with catastrophes was considered in Krinik and Mortensen 
(2007). 

Models based on stochastic processes in the presence of catastrophes have been recently 
exploited also in another field of mathematical biology, with special reference to the description 
of the interaction between a myosin head and an actin filament. Precisely, in Buonocore et 
al. (2009) a continuous-time Markov chain-based model and its continuous approximation have 
been proposed and studied starting from the M/M/l/K-catastrophe paradigm. 

Along the research line traced by the above mentioned articles, in this paper we consider a 
different birth-death stochastic model subject to catastrophes. We recall that a bilateral birth- 
death process characterized by constant birth and death rates, and subject to jumps to state 
induced by catastrophes, has been studied recently in Di Crescenzo and Nastro (2004) and in 
Switkes (2004). Here we aim to study an extension of such a stochastic model, which performs 
a continuous-time random walk on the integers and, in addition, is subject to the occurrence of 
catastrophes followed by exponentially-distributed repair times. The state-space of this model 
thus consists in the whole set of integers complemented by a spurious state, say F, which is 
occupied by the system during the repair times. At the end of any repair the system occupies 
the state zero and starts afresh. It should be mentioned that the idea of an M/M/1 queueing 
system subject to catastrophes and followed by random repairs has been recently considered in 
Krishna Kumar et al. (2007). 

It is worth pointing out that the basic stochastic process of our model, namely the bilateral 
birth-death process with constant birth and death rates, can be viewed for instance as describing 
the state of a double-ended queue in a taxi-passenger system. Double-ended queues have been 
studied e.g. in Conolly et al. (2002). Here, denoting by {N{t); t > 0} the state of the considered 
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system, N{t) = n > means that at time t there are n customers in the system, A^(t) = — n < 
means that at time t there are n taxis in the system, N{t) = means that at time t the system 
is empty, and N{t) = F means that at time t the system is in a repair period. 

The transient analysis of queueing systems is not always easy to be performed, and it leads 
very often to non-manageable formulas. In many cases one is forced to study the system under 
heavy-traffic conditions, in order to obtain at least diffusion or continuous approximations. Also 
systems subject to catastrophes can be approximated in such a way, the approximating process 
being usually a jump-diffusion process. This is the case of the contribution given in Di Crescenzo 
et al. (2003), where an M/M/1 queueing system in the presence of catastrophes is studied and 
is approximated by a Wiener process subject to randomly occurring jumps. A similar result will 
be obtained hereafter for the system under investigation. 

This is the plan of the paper. In Section 2 we study the preannounced double-ended queue 
with catastrophes and repairs. By exploting its connection with the system in absence of catas- 
trophes wc express the state probabilities of the queue in terms of integrals of Bcssel functions. 
The mean and the variance are also investigated, together with the steady state of the system. 
Section 3 is devoted to develop the approximating procedure that leads to a jump-diffusion pro- 
cess of the Wiener type. We thus analyze the approximating process, and obtain its transient 
probability density, mean, variance, and steady-state density. In conclusion, in Section 4 we 
discuss the goodness of the continuous approximation. 

2 Discrete model with catastrophes and repairs 

We consider a system performing a one-dimensional random walk in continuous time on the 
whole set of integers. Let {N(t); t > 0} be the position of the system at time t. The system 
moves along the real axis such that A is the rate of moving to the right and /x is the rate of 
moving to the left (see Conolly, 1971). This random walk process {N{t); i > 0} is also said 
to be a bilateral birth-death process, with birth rate A and death rate in any small interval 
{t,t + At), At > 0, a birth occurs with probability XAt + o{At); a death occurs with probability 
IJ,At+o{At). It is clear that in this interval neither a birth nor death takes place with probability 
1 - {\ + n)At + o{At). 

Apart from birth and death processes, catastrophes occur at the system according to a 
Poisson process with rate i.e., the catastrophe occurs at the system in the small interval 
{t, t+At) with probability iyAt+o{At). Whenever a catastrophe occurs at the system, the system 
goes into the failure state F, i.e., the bilateral birth-death process is subject to catastrophic 
failure. The repair times of the failed system are i.i.d. according to an exponential distribution 
with mean 7]^^. After a repair of the system is completed, the system immediately returns to 
state zero and begins to work again ("on" state). 

As mentioned in Section 1, N{t) can be seen as the state of a double-ended queue subject 
to disasters and repairs. We model the system as a continuous-time Markov process with state- 
space S = {F}UZ = {F,0,±l,±2,...}. Let P„(t) = P[N{t) = n], n = 0,±1,±2,..., be the 
state probability that the system is in state n at time t when the system is working ( "on" state) 
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and q{t) be the probabihty that the system at time t is under repair ("failure" state). 

For the system under investigation, the Chapman-Kolmogorov forward differential-difference 
equations for the system state probabilities Pn{t), n = 0, ±1, ±2, and failure state probability 
q{t) can be written as 

dq{t) 



dt 
dPo(t) 
dt 

dPnit) 



r]q{t)+u[l-q{t)], (2.1) 
-(A + /i + i^) Po{t) + A P_i(t) + fiPi{t) +iiq{t), (2.2) 
-{X + fi + u)Pn{t)+XPn-l{t) + flPn+l{t), n = ±1, ±2, ±3, ... . (2.3) 



dt 

Without loss of generality, assume that initial position of the system at time t = is zero, i.e., 

q{0)=0 and P„(0) = P if ™ = (2.4) 

L otherwise. 

2.1 Transient distribution and relation with the process in absence of catas- 
trophes 

Let us now introduce the so-called randomized random walk {N{t); t > 0} with parameters A 
and /i. It is a bilateral birth-death process over Z characterized by birth rate A and death rate 
^, and describes the state of the system under consideration in the absence of catastrophes. We 
recall that the state probabilities of process N{t) conditional on A^(0) = are given by (see, for 
instance, Baccelli and Massey 1989 or Conolly 1971) 

P^{t) = e-(^+/^)*/3" /„ (at) , t > 0, n e Z, (2.5) 

where 

a = 2yV, /3 = ^ (2.6) 

and /«(•) is the modified Bessel function of the first kind of order n. We remark that the 
distribution given in (|2.5p is often named Skellam distribution. 

Hereafter we express the state probabilities Pnit) in terms of suitable quantities concerning 
the bilateral birth-death process in absence of catastrophes N{t). Indeed, from probabilistic 
arguments it is not hard to see that the probabilities Pn{t) can be expressed in terms of (|2.5p 
as follows: 

Pn{t) = e-"' Pn{t) +11 f q{T) e-^ P„(t - r) dr, t > 0, n G Z. (2.7) 

JO 

The first term in the right-hand-side of (j2.7p is the probability that the system goes from to n 
without occurrence of catastrophes in (0, t), whereas the second term represents the probability 
that at time r G (0, t) the system enters state due to a repair and then goes from to n 
without occurrence of catastrophes in (r, t) . 



Theorem 2.1 For all t > we have 



lit) 



I _ g-(^+'^)*' 



4 



-ut 



-''^-(^+'^)(*-^)4(a(t-T))dr 



n G 



(2.9) 



where a and j3 are defined in (E 



Proof. Making use of Eqs. (12. ip and recalling the first of (12. 4p we obtain (12. Sp . Furthermore, 
substituting ([23]) and ([M]) in Eq. ([12]) the transient distribution ([22]) follows. □ 

Some plots of probabilities (12. 9p are shown in Figure 1. The expressions given in Theorem 
12.11 completely determine all state transient probabilities of the system and the failure state 
probability q{t) for the bilateral birth-death process N{t) with catastrophes and repairs. 

Let us now denote by T^ n the first-passage time from state to state n {n = ±1, ±2, ±3, . . .) 
of process N(t), and by 5o,n(i) the corresponding probability density function. It is well known 
that (see Abate et al. 1991, for instance), for all t > 



go At) = -fPnit) = ^-fe-^^^^^'p^In{at), n = ±1, ±2, ±3, . . . 



(2.10) 



Due to the Markov property, the following renewal-type equation holds for the system under 
consideration, for t > 0: 

r-t 



Pn{t) 



Po(^^)e-"(*-")5o,n(i-n) du, 



n = ±1,±2,±3, ... 



(2.11) 



Hence, Eq. ()2.1ip expresses that any sample-path of N[t) going from to state n, at an eventual 
time u € (0,t) visits the state and then goes from to n without occurrence of catastrophes 
in (n, t). Hereafter we provide an alternative expression of the transition probabilities. 

Proposition 2.1 For all t > the following relation holds: 



Pn{t) = |n|/3" / e 

'0 



t — U 



n = ±l,±2,±3,... 





0.0 0.2 0.4 0.6 OS 1.0 1.2 1.4 



Figure 1: Plots of Pn{t) for n = 0,1,2 (from top to bottom near the origin), where fi = 2, 
V = 0.1, r/ = 1, and (a) A = 2, (b) A = 8. 
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Proof. It immediately follows by making use of (j2.10p in (j2.1ip . □ 

In the following, for any function f{t), t > 0, we denote by f*{z) = e~^^ f{t) dt, z > 0, 
its Laplace transform. 

By adopting customary methodologies, such as those described in Jain et al. (2007) or 
Parthasarathy and Lenin (2004), we can evaluate the Laplace transform of the transition proba- 
bilities Pn{t)- Indeed, by taking Laplace transforms in the system of equations (|2.1|) - (j2.3|) . after 
some calculations we obtain: 

P^{z) = ^ + , "^"^ (2.12) 

^J{z + A + ^ + zy)2 - 4A/X z{z + r} + v)^J{z + A + /i + i/)2 - 4A^ 



and 



where 



P*i,) = \ P5{^)[M^T^ n = 1,2,3,... 

1 Po*(z)[^i(z)]-, n = -l,-2,-3,... 



I t ^ , t ^ z + \ + ii + u±^{z + \ + ii + vf I ( I ( ^ 

i)l{z),llj2{z) = , ll)l{z) > '>p2[z) 

are the roots of quadratic equation fix^^z) — (z + A + /i + v)x{z) + A = 0. Expressions (|2.12p and 
(j2.13p will be used in Section 4 in order to discuss the goodness of the approximation of N{t) 
by a suitable jump-diffusion process X{t). 

2.2 Steady-state distribution 

In this section, we shall investigate the structure of the steady-state probabilities and the failure 
state probability of process N{t): 

Tin ■■= lim P„(t), q := lim q{t). 



Theorem 2.2 The steady-state probabilities of the system state {7r„ : n = 0, ±1, ±2, . . .} and 
the failure probability q of the bilateral birth-death process with catastrophes are obtained as, for 
u > and r] > 0, 

(2.14) 
(2.15) 



rj-\- ly 

TTo = (1 - ( 



TTr, 



2ii 



n = 1,2, 



(2.16) 



and 



IT- 



{X + n + iy)- {X + + - AXfi 



2A 



TTo n = l,2, .... 



(2.17) 
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Proof. The steady-state probabilities of the system can be obtained by letting t — )• +00 in 
the system of equations (|2.1l) - (l2.3p . The resulting balance equations yield a system of recur- 
rence relations of degree 2 with known coefficients, whose solution can be obtained by standard 
methods. Eqs. (fmi) - (f2l7]l then follow. □ 

2.3 Moments 

Let us now focus on the mean of the bilateral birth-death process with catastrophes and repairs, 
and set 

mjv(t) = E[Nit) ■ l{Nit)^F} 1 A^(0) = 0], t>0. 
Then, from (j2.7p we have 



TUNit) = e-^^mit) +r] [ q{T) e'" ^^'^'^ mj^{t - t) dr, 

Jo 



t > 0, 

'0 



where m^{t) := E[N{t) \ N{0) = 0] = (A — fi)t is the mean of the randomized random walk with 
initial state 0. Hence, recalling ()2.8p after some calculations we obtain 



mN{t) = ^— ^ 1 - e-^* + -e-''\l - e"^*) . (2.18) 
[T] + v)v IT J 

Due to (I2.18p . if A > the following holds: 

- if r/ > z/ then miy{t) is increasing for all t > 0, 

- a rj < u then m (t) is increasing for t < t and is decreasing for t > t, where 

- 1 , 
t := - log' 



r] V — 7] 

Clearly, if A < the monotonicity of mf^{t) is reversed. Some plots of mf^{t) are shown in 
Figure 2. 

Let us now introduce the variance of the bilateral birth-death process with catastrophes and 
repairs: 

VN{t) = Var[N{t) ■ l{N{t)^F} I N{Q) = 0], t > 0. 
By adopting a similar procedure, making use of (j2.7p . ()2.18p and recalling that Var[N{t) \ N{0) = 
0] = (A -|- /i)t, it is not hard to obtain VAr(i). Some plots of the variance of N{t) are given in 
Figure 3. 

We now discuss the asymptotic mean and variance of process N{t). From (|2.18p it is not 
hard to see that, for v > 0, 

lim m^it) = = (1 - g) (2.19) 

We point out that the presence of catastrophes and repairs may be seen as a regulatory effect 
to the system. Indeed, if = the system does not admit a steady-state distribution, whereas 
if > then the right-hand-side of (I2.19P is finite. Finally, straightforward calculations lead to 
the following asymptotic variance, for w > 0: 

lim V^{t) = + ,^^^(2. + ,). (2.20) 
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3 Jump-diffusion approximation 



Let us now construct a continuous approximation starting from the process N{t). We shall make 
use of a scaling similar to one succesfully used elsewhere in queueing theory contexts (see Giorno 
et al. 1986, 1987; Di Crescenzo and Nobile 1995). We recall that a popular paper dealing with 
Brownian motion with jumps that approximates a standard single-server queue with random 
vacations is due to Kella and Whitt (1990), whereas Kimura (2004) proposed an useful review 
of some basic concepts and issues in diffusion modeling, and a bibliographical guide to diffusion 
models for queues that arc typically found in computer/communication systems. 

Let us consider the Markov process {N^{t); t > 0}, having state space {F}U{. . . , — e, 0, e, . . .}, 
with £ > 0. For all t > 0, the transitions of N^(t) are governed by the following rates: 



lim -^PiN^it + At) 
Atio At ^ ' 



(n + l)£|Ar«(t) 



n G Z, 
n G Z, 



where A, /i and a are positive parameters. Moreover, the process \N'^{t)\ t > 0} is still subject 
to catastrophes arriving with rate v. Repairs arc still exponentially distributed with rate and 
the system enters state just after a repair. Hence, the repair and failure time distributions are 
not affected by the scaling procedure introduced above, so that 

A40 i^^^^"^* + At) = F I K(t) = ne) = v, n G Z, 



lim ^P{N^{t^ At) 



where rj and u are the same parameters of N{t). Essentially, the integers of the state-space of 
N{t) have been rescaled by a factor e, and the following substitutions have been performed: 



, _ A a' 

^-7 + 2i^' 



(3.1) 



rriNif) 





(b) 



10 20 30 40 50 60 70 



Figure 2: Plots of mN{t) for A — = 1, 77 = 0.25, 0.5, 1, 2 (from bottom to top), and (a) = 1, 
(b) ly = 0.1. 
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For t > the transient behaviour of N^(t) is described by the following probabilities: 

P^{t) = P{N^{t) =ne}, nGZ. (3.2) 
As before, assuming that P{N^{0) = 0} = 1, we have Pq{0) = 1. Moreover, we note that 

^P„"(t) = l-g(t), 

nGZ 

where q{t) is the probability that N^{t) is in state F, and is given in (j2.8p . 

Let us now consider a continuous-time Markov process {X{t); t > 0} having state space 
{F} U (— oo,cxd). This is a jump-diffusion process that describes the motion of a particle that 
starts from the origin at time 0, then behaves as a Wiener process with drift A — /x and infinitesi- 
mal variance cj^, until it jumps to state F with rate i^. Subsequently it jumps again to the initial 
state 0, with rate rj, and diffuses afresh. The probability law of X(t) possesses an absolutely 
continuous component 

f{x,t\0) = \iml-P{x < X{t) < x + h\X{0) = 0}, x G M, t>0 

and a discrete component P{X(t) = F \ X{0) = 0}. The latter is easily seen to be identical to 
the probability q{t), obtained in (j2.8p . 

The continuous approximation relies on the convergence of N^{t) to process X{t) as e tends 
to 0. Probabilities ()3.2p are near to f{n£,t\0)e when e is close to zero, and the correspon- 
dence between such quantities becomes increasingly better as e approaches zero. We thus follow 
a customary procedure based on the substitution of Pn{t) with /(ne,t|0)e in the Chapman- 
Kolmogorov forward differential-difference equations of N^(t). After setting x = ne and ex- 
panding / as Taylor series, we pass to the limit as e J, 0. By taking positions (|3.ip into account 
and assuming that At is proportional to e^, after some calculations we obtain: 

^/(x,t|0) = -z./(x,t|0)-(A-/I)^/(x,t|0) + y ^/(x,t|0), (3.3) 
to be solved with condition 

/oo 
f{x,t\0)dx = l-q{t), t>0. (3.4) 
-oo 
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Equation (|3.3p is the Fokker-Planck equation for a Wiener process perturbed by a leaky term 
—vf{x,t I 0). Note that the initial condition (0) = 1 now becomes 

lim/(x,t |0) = (3.5) 

where 5{x) is the Dirac delta- function. 

The above sketched procedure indicates that the process N^{t) converges in some sense to 
the process X{t) as e tends to 0. Hence, statistical indexes of X{t) are suitable to determine 
quantities useful to describe N{t) when e is very small. Indeed, under position (|3.ip we have 

P{N{t) < n I N{Q) = 0} ~ P{X{t) <ne\ X(0) = 0}; 

recalling that x = ne, this approximation is expected to improve as e goes to zero and as \n\ 
grows larger. 

3.1 Transient and steady-state distributions 

In order to analyse the process X[t) hereafter we obtain the function /*(x,z|0), that is the 
Laplace transform of density | 0). Recalling (|3.5p . from (j3.3p we then have 

(z + z.)r(x,z[o) = -(A-/i)^r(x,z|o) + y^r(x,z|o), x/o, (3.6) 



whereas condition (13.41) becomes 



r(x,z|0)dx= /^"^ (3.7) 



Theorem 3.1 For z > we have 

1 {z + i^){z + r]) 



r(x,z|0) 



A_M A/(A-/2)2 + 2a2(z + z.) 

X exp { ^^-2^ X - ^ ^|x|} (xGM). (3.8) 

Proof. By continuity of /*(x, 2: | 0) in x = the solution of ()3.6p takes the form 

r(x,.|0) = |^^^^^ (3.9) 
^ ' ^ Me"'!^ ifx<0 

where 



A-/I ± JiX-fiy + 2a^z + i^) 

Wi,W2 = ^2 {W2<0<wi) (3.10) 



a 

are solutions of equation 

.2 ...2 



- 2(A - /2) w - 2(z + i^) = 0. 
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Making use of (j3.9p in (j3.7p we thus obtain 



(A-/I)2 + 2^72(z + r/) 



Hence, Eq. dM]) finally follows from (IXTIl and (pTTT]) . 



2; (z + 77 + z^) 



(3.11) 



□ 



We are now able to obtain an expression of | 0), that is the analogous of (j2.7p for the 

discrete system. 



Theorem 3.2 For all x gR. and t > we have 



/(x,t|0) = e-'^*/(x,t|0)+77 / g(r)e 

Jo 



/(x,t-r|0) dr, 



(3.12) 



where f{x,t \ 0) is the transition density of a Wiener process {X{t); t > 0} having drift X — Jl 
and infinitesimal variance . 

Proof. First of all we note that Eq. (j3.8p can be expressed as 

rju 1 rju 1 



f*{x,z\^) 



, - ,^ - \ A -/I 

1 H exp < — ^ X > 

rj + u z rj + uz + r] + vj I cr^ I 

1 I J{X-f,Y + 2a^{z + v) 



{\-fif + 2<j^{z + v) 



: exp 



0"^ 



(x G R). (3.13) 



Let us now consider a Wiener process having drift y (A — /x)^ + 20"^^ and infinitesimal variance 
(T^; its transition density fw{x,t \ 0) has Laplace transform 



fw{x,z\Q) 



(A-/I)2 + 2a2(z + z.) 



'(A-/2)2 + 2a2z. 
exp < T\ X 



{\-ri? + 2a^{z + u) 



so that 



exp <^ V (A - A*)^ + + ^) 



= exp 

'(A-/2)2 + 2a2(z + i.) 

Making use of Eq. ([314]) . thus Eq. (f3J3]l yields 

rjv 1 r]v 1 



(A - /2)2 + 2ct2i/ 



xW;-(x,z|0) (3.14) 



/*(x,z|0) =e^2x 1 + 



r] + u z rj + uz + rj + u 



fl{x,z\0) (xGM), 



(3.15) 



where we have set 



Z2 



J{\-fiY + 2a^v 



Hence, by taking the inverse Laplace transform in (|3.15p . we come to 



f{x,t\d) = e'^'^ lf^{x,t\d) + 



r] V 
V + ^ Jo 
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1 - e 



fwix,T\ 0) dr 



Finally, recalling (|2.8p and making use of identity 



oZ2 X 



fw{x,t\ 0) 



-vt 



/(^,t|o), 



Eq. (IXT2]) follows. 



□ 



Some plots of density /(x, f | 0) obtained from (|3.12|) are given in Figure 4. 
Hereafter we obtain another expression of /(x, t \ 0) which is the analogous of (j2.1ip under 
the continuous approximation. It involves the first-passage-time density of the Wiener process 
with drift A — /I and infinitesimal variance cr^, given by: 



g[x,t\xQ) = /(x,t|xo), x/xo, t > 0. 



(3.16) 



Theorem 3.3 For all x eR \ {0} and t > we have 

fix, t\0)= f /(O, r I 0) e-^ (*-^) g{x, t - r | 0) dr. 
Proof. We recall that for all x G M and t > the following relation holds: 



/(x,t|0) =exp<; ^^^^xi /(0,t|x). 



(3.17) 



Hence, for x € M and t > Eq. (I3.12p becomes 



/(x,t|0) = e~'^'/(x,t|0)+r? / g(r) e 

JQ 



- it-r) J X i /(O, t - r I x) dr. 



(3.18) 



We recall that for the Wiener process X{t) one has: 



/(0,t|x)= / 5(0,r|x)/(0,t-T|0)dr, x 7^ 0, t > 0. 
^0 



(3.19) 



f(x,t\0) 


0.35 




0.30 




0.25 




0.2(Y 




o/l5/ 




/o/io 




//o.pi 





Figure 4: Plots of density (|3.12p . for A - /J = 2, r/ = 1, cr^ = 1, t = 1, and v = \, 0.5, 0.1 (from 
left to right). The probability mass is respectively 0.5677, 0.7410, 0.9394, due to ([37 
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Moreover, making use of (j3.16p and (j3.17p . for x ^ and t > we have 



5(0, 1 1 = M exp I -^^^^ X \ fix, t \ 0) 



exp<( - ^''"^ ^ x \ g{x,t\0). (3.20) 



Substituting (13.19P in the last term of (I3.18p . employing Fubini's theorem and making use again 
of (j3.12p we finally obtain: 

fix,t\0) = e-''^f{x,t\0)-e-''* exp | ^-^^^^ x | ^ g{0,u\x) fiO,t - u\0) du 

+ exp I ^ x| ^ g(0,tt|x)e~'^"/(0,t-n|0)dn, x / 0. (3.21) 

The first and the second term in the right-hand-side of (j3.2ip vanish due to Eqs. (j3.17p and 
(|3.19p . Hence, the thesis follows making use of ()3.20p in the last term of (j3.2ip . □ 

Aiming to analyse the asymptotic behaviour of process X{t) we now introduce the steady- 
state density 

W{x)= lim f{x,t\0), xGR. (3.22) 

t— s>+oo 

Hereafter we obtain that W{x) is a bilateral asimmetric exponential density. 
Theorem 3.4 The steady-state density i3. 22\) is given by 



Wix) = ^ = exp J X - /x)^ + 2a2z. ^ ^ 

^ + ^^(A-/2)2 + 2a2^ i ^ 

Proof. Being W{x) = \im\z f*{x, z\0)], from ([HS]) we have ([323]). □ 

2:4.0 



We remark that, due to (j3.4p and (j2.14p . we have 

/oo 
W{x) dx = 1 — q 
-00 



7] + U 



3.2 Mean and variance 

Let us now define the mean 

mx{t) := E[X{t) ■ l{x{t)^F} I ^(0) =0], t>0. 
Due to (|3J2]1 and recalling that E[X{t) \ X{0) = 0] = (A - Jl)t, for all t > we have 

mx (t) = e-"' (A -J2)t + r] f q{T) e"'^ (A - /I) (t - r) dr. 







13 



mx(t) Vxit) 
20.0 r 




Figure 5: Plots of mean and variance (j3.24p and (|3.25p for the same values of the parameters as 
in Figure 4 (from bottom to top). 



Use of ()2.8p and straightforward calculations thus yield: 

(A-/2)r? r ^ 1/2 



mxit) = , ^ ' 1 - e-'^' + -^e-^'il - e-^') , t > 0. (3.24) 

The strict correspondence between the means (I2.18P and (I3.24p is evident. 

Making use of (I3.12p one can express the moments of X{t) in terms of the moments of the 
Wiener process. In particular, for any t > we have: 

Vx{t) := Var[X{t) ■ I ^(0) = 0] 



1 - e"^* + -2 e-*^* (1 - e-'?*) 



+ r/^ + 2riu{ri + u){u^ - jf)te~''' - e'^'"' (i/^ - rj^ - i/^e-"*)^ |. (3.25) 

Notice that for the densities shown (from left to right) in Figure 4 we have mx(l) = 0.865, 1.305, 
1.834 and Vx(l) = 1.283, 1.469, 1.198 (remarking that in this case Vx(l) is not monotonic in 
v). The behaviour of mean and variance as t varies is shown in Figure 5. 

Finally, from (j3.24p and (j3.25p we immediately obtain the asymptotic limits of the mean and 
variance of the approximating process: 

lim mx(i) = |^^, lim Vx{t) = -^V + + (^.26) 



4 Remarks on the continuous approximation 

We conclude the paper by giving some comments concerning the goodness of the continuous 
approximation. Recall that, as pointed out in Section 3, the approximation is expected to 
improve as e goes to zero and as |n| grows larger. 
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First of all we point out that the means of the discrete process and of the continuous ap- 
proximating process are in strict agreement. Indeed, by comparing (j2.18p and (I3.24p . under 
substitutions (13. ip we immediately have 

-X{t) XiO) 



E[N{t)\N{0) =0]=E 







for all e > and t > 0. 

In order to ascertain the approximation between the transient distributions of N{t) and X{t), 
let us now denote by P*{z;e), n E Z, the Laplace transform of Pn{t) where the substitutions 
p.ip are performed. Hence, making use of (j2.15p and (|3.8p we obtain 

lim5L(£i£)= ^r(o..|o). 

Moreover, by noting that under positions (j3.ip 



z + \ + ^ + v - J{z + X + fi + uy -4Xfi ^ 
lim = 1, 

from (I2.13P for x ^ we have 

j e^^-riO,z\0), x>0\ 

hm — = < > = f (x, z \ U), 

eio e y e"'i^/*(0,z|0), X < J 

where f*{x, -z | 0) is given in (|3.8p and where wi and W2 are defined in (j3.10p . Hence, recalling 
the meaning of P*[z;£) and denoting by Pn{t;£) the probabilities ()2.9p when the substitutions 
p.ip are performed, we have: 

lim ^"^^'^^'"="/^ =/(:r,t|0). 

This implies that under the approximating procedure the state probabilities of the discrete 
process are close to the density f{x,t \ 0) of process X{t). In addition, this is confirmed by the 
data of Table 1 showing that the agreement of the two steady-state distributions improves as e 
decreases and as \n\ increases. Furthermore, we remark that performing the substitutions 

£h = £h, Xh = Xh, 'j2h = 'flh, af^ = a'^h, (h > 0) (4.1) 

recalling (j3.ip we have 



Eh 2e^ £h 2£l 



Hence, carrying out substitutions (14. ip in the steady-state density (I3.23P we obtain Wh{n£h)£h = 
W{n£)£ for any h > 0. This shows that the approximation procedure is effective when the rates 
A and fi are suitably varying. 

In conclusion we point out that the agreement between the discrete model and the approx- 
imating jump-diffusion process is also confirmed by the comparison of their asymptotic means 
and variances. Indeed, the correspondence of limits ()3.26p with those obtained in Eqs. (j2.19p 
and (|2.20p for the discrete model is immediate under positions ()3.ip . 
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£ 


= 0.1 










e 


= 0.05 










e 


= 0.01 






n 






W{ne) 


A(n) 


TT 




W{ne) 


A(n) 


TT 




W{ne) 


A(n) 


-6 





03621 





03668 





01305 





04073 





04102 





00721 





04480 





04487 





00155 


-5 





03757 





03807 





01353 





04149 





04180 





00733 





04496 





04503 





00156 


~4 





03897 





03952 





01401 





04227 





04258 





00745 





04513 





04520 





00156 


-3 





04044 





04102 





01448 





04306 





04339 





00757 





04530 





04537 





00157 


-2 





04196 





04258 





01497 





04387 





04420 





00769 





04547 





04554 





00157 


-1 





04353 





04420 





01544 





04469 





04503 





00781 





04564 





04571 





00158 








04516 





04588 





01593 





04552 





04588 





00793 





04581 





04588 





00158 


1 





04260 





04323 





01472 





04420 





04454 





00763 





04554 





04561 





00157 


2 





04019 





04073 





01351 





04292 





04323 





00732 





04527 





04534 





00156 


3 





03791 





03838 





01231 





04167 





04196 





00702 





04500 





04507 





00154 


4 





03576 





03616 





01111 





04046 





04073 





00672 





04473 





04480 





00153 


5 





03373 





03407 





00990 





03929 





03954 





00641 





04447 





04454 





00152 


6 





03182 





03210 





00870 





03814 





03838 





00611 





04421 





04427 





00151 



Table 1: For rj = 0.25, z/=l, A = l, /U = 2 and a = 3 the quantities 7r„(e)/e and W{n£) are listed 
together with the relative difference A(n) := [W{ne)e — 7r„(e)]/7r„(e) for various integer values 
of n and various choices of e. We have (A,/Lt) = (460,470) when e = 0.1, (A,/x) = (1820, 1840) 
when e = 0.05, and (A,/x) = (45100,45200) when e = 0.01. 
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